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Abstract 

The interaction of coherent magnetization rotation with a system of two-level 
impurities is studied. Two different, but not contradictory mechanisms, the 
'slow-relaxing ion' and the 'fast-relaxing ion' are utilized to derive a system 
of integro-differential equations for the magnetization. In the case that the 
impurity relaxation rate is much greater than the magnetization precession 
frequency, these equations can be written in the form of the Landau-Lifshitz 
equation with damping. Thus the damping parameter can be directly calcu- 
lated from these microscopic impurity relaxation processes. 



1 



Typeset using REVTfjX 



I. INTRODUCTION 



The problem of magnetization reversal in a single domain ferromagnetic grain is closely 
connected to the problem of relaxation. The channels of energy transfer from the coherent 
magnetization rotation to the thermal bath can be defined by interactions with itinerant 
electrons, spin waves (magnons), lattice vibrations (phonons), and crystal imperfections 
(impurities, defects). The theories of magnon damping in itinerant ferromagnets give 
relaxation rates oc k 2 , where k is the wave vector. Therefore this channel of energy loss is 
expected to be very small for coherent magnetization rotation (k = 0). Very large mag- 
netization motions are accompanied by nonlinear magnon excitations |||4[]. However it has 
been shown that this mechanism is ineffective in extremely small grains that undergo uniform 
magnetization reversal ||. A mechanism for direct coupling to the lattice via magnetoelastic 
interaction has been proposed ||. 

Here we focus on the impurity relaxation processes. We believe that in fine grains the 
role of imperfections is expected to be much larger than in bulk material. In addition to 
paramagnetic ions (such as Fe 2+ , Fe 4+ , Mn 3+ , Cr 2+ , Cr 4+ etc.), the surface magnetic atoms 
are expected to behave as impurities. 

It is well known that impurity magnetic ions of transition and rare-earth elements with 
strong spin-orbital coupling give rise to ferromagnetic resonance line broadening ||[7] . Two 
different (but not contradictory) mechanisms of relaxation have been proposed. In the, so- 



called, 'slow-relaxing ion' mechanism |qjl]]--[l5| small magnetization oscillations modulate 
the impurity levels in the vicinity of thermal equilibrium. This energy modulation is ab- 
sorbed by the lattice at an extremely fast rate. The 'fast-relaxing ion' mechanism [|9|, p~Q|JT~2 



is effective when the transverse relaxation rate of the impurity is so fast that the magnetiza- 
tion oscillations can excite direct transitions between impurity levels. We present a theory of 
magnetization relaxation of impurities for large magnetization oscillations, including mag- 
netization reversal. 
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II. MODEL 



The fundamental magnetization dynamics is described by gyromagnetic precession: 

— = -7m X H eff , (1) 

where, for coherent rotation, H e g- = —dU/dWl is the effective magnetic field defined by the 
energy density U, M = M s m, (|m(£)| = 1) and 7 is the gyromagnetic ratio. We shall 
consider the impurities as independent two-level systems (Fig.l). Let us write the energy 
density in the form U = U m + Ui mp , where W m is the magnetic energy density (including only 
anisotropy and Zeeman energies). The impurities energy density can be written in the form: 

W imp = (h/V) ^0,j n j + W ox, 
3 

W ex = y Et w i • m " V2) + e, • m (c+) + e* • m ( Cj )}. (2) 
3 

where j indicates a sum over all impurities. 

As illustrated in Fig.l, V is the grain volume, Q j is the splitting (due to crystal fields, 
exchange and spin-orbit interactions) of the j-th impurity. Wj, ej (and its complex conjugate 
e*) characterize the anisotropic exchange interaction between each impurity (as an effective 
spin 1/2) and host nearest neighbor spins (derived in the Appendix). In Eq.(|2]), rij = (CjCj 
is the upper level impurity population, operators Cj and Cj describe transitions between the 
levels and (...) denotes quantum mechanical expectation values. The contribution of the 
impurity interactions to the effective field is: 

rfl/l h 

^ = " MV £ - V2) + ^' + ^ • (3) 

s 3 



III. EQUATIONS OF MOTION 

The Eqs.(l),(3) describe the dynamic magnetization including coupling to the impurities. 
The interaction of j-th impurity with the thermal bath will be considered phenomenologi- 
cally, introducing the longitudinal and transverse T± t j relaxation rates which can depend 
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on the impurity splitting Qj, temperature T and external static magnetic field Ho. We also 
assume that the characteristic frequencies of magnetization rotation are much smaller than 
Qj. Thus the kinetics of the j-the impurity population can be defined by the following 
equation: 

J t n i = - r |],i(^i)K' - M^j))- ( 4 ) 

It is assumed that at each instant of time the population rij relaxes to the equilibrium value 
riT^j) = [exp^hQ/ksT) + l] -1 corresponding to the dynamic splitting Qj = Qqj + SQj(t), 
where 5Qj(t) = Wj ■ m = Wj tX m x (t) +Wj ty m y (t) +Wj yZ m z {t). Note that the rate Yy{Vtj) also 
depends on this dynamic splitting. The kinetics of "transverse" impurity spin components 
is given by: 

(cj) = -iVLj (cj) - iej • m [(1 - 2 nj )]. (5) 

A. General solution 

The general solution of Eq.(|4]) is: 

t 

5n 3 {t) = J exp[$ y (tO - $ y (jt)] ryinjfaMnrpjfa)] - n T (^ j)}^i, (6) 

— oo 

where 

t 

$ y (t) = J r y [^(t')]< 5n 3 (t) = nj (t) - nr(noj). (7) 

— oo 

Analogously, the solution of the Eq.(||) can be written as: 

t 

( Cj (t)) = -i J exp[$ ±J (*i) - $±j(t)] ej ■ m(ti) [1 - 2n T (n ,j) - 26n j (t 1 )]dt 1 (8) 

— oo 

with 

t 

*±j(t) = J {rx^PMOl + i^(t')}dt'. (9) 

— oo 

By substitution of Eqs. (§),(§D into the effective field ([$]) and the latter into the Landau- 
Lifshitz equation ([[]), the system of integro-differential equations for the magnetization m 
is obtained. 
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B. Slow relaxation 



In the case of small modulation h\8Vlj{t)\/k-QT <C 1 we can write the following expan- 
sions: n T [Q j + 5Qj(t)} = nr(fioj) + [dnri^oj) / dft 0t j]8Qj(t) + ... and ry[O ,.j + 5Qj(t)] = 
T\\j(floj) + [dT\\j(Q j) /dQo t j]5Qj(t) + ... Utilizing only the first expansion terms, Eq.(Q) 
becomes 

8nf\t) = TyiQoj) dn g^°; j) } exphrnjCOoj)^ - h)) Wj • mfa) dt x . (10) 

J -oo 

For the case of short memory (longitudinal relaxation T\\j(Qqj) is much faster than the 
characteristic frequency of magnetization rotation) we can introduce a variable r = t — t x 
and make a moment expansion in powers of r. The result has the form 



1 dm(t) 

m{t) ~r^\~) dt + 



(ii) 



This first order expansion to the kinetic equation (4) for the impurity relaxation yields an 
additional term to the effective magnetic field @: 



<9n<r(^o,j) 1 dm(£) 



-}• (12) 



The first and the second terms in Eq.(|i"2"D can be included into the effective anisotropy fields. 



The term containing dm(t)/dt in ( |i2|) gives the magnetization relaxation 



C. Fast relaxation 

The kinetics of the transverse impurity components @ in the simplest approximation 
can be written as: 

t 

( Cj (t)) ~ - 2n T (Sl 0ij )] / exp[(r ±J + - *)] e,- ■ m(^) dt x . (13) 

— oo 

For a short memory the moment expansion gives: 



<<=,(«» - - i i - 2 "rS oj> « 



m(t) 



dm(t) 



(14) 



This solution (and its complex conjugate (c+(t)y = (cj(t)}*) may directly be added to the 
effective field (3). 



IV. EXAMPLE 

Consider for simplicity Qqj = and the case of slow relaxation where Wj = (w,0, 0) 
and ej = 0. Let us write the magnetic energy in the form U m = K u (l — m z 2 ) — HoM s m z , 
which assumes an applied field H parallel to a uniaxial anisotropy axis. In this case the 
combination of Eqs.(l),(3) and (12) may be written as: 

dm x 

— — = -{m z + h )m y , 

CLT 

-j-^ = {m z + ho)m x + a 1 (dm x /dr)m z , 
dr 

—± = - ai (dm x /dT)m y , (15) 
where r = ^H^t, h = H /H K , H K = 2K U /M S and 

ai - c imp ^ J ^ (nflo/feT) + 1]2 > c imp - M ^ . U bj 

iVi m p is the number of impurities in the grain and z- imp is the average number of magnetic 
neighbors for one impurity. 

Fig. 2 shows a typical reversal of magnetization described by the Eqs.([[|). It is interesting 
to note that the Landau-Lifshitz equation with phenomenological damping 1Z = —a 7 m x 
(m x H e ff) gives virtually the identical evolution (within 1% of accuracy) of m z (r) if a = 
ai/2. This correspondence between a and oti is exact for small magnetization motions. 

This work was partly supported by matching funds from the Center for Magnetic Record- 
ing Research at the University of California - San Diego and CMRR incorporated sponsor 
accounts. 
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APPENDIX A: DERIVATION OF EQ.(2) 



We assume that the exchange interaction between the impurity and the neighboring host 
atoms can be written in terms of an effective spin (s = 1/2) of the two-state impurity ion. 
In this case the Hamiltonian may be written in a general form: 

We* = E E ^M R i' r^Sa^SaiRj + T„). (Al) 

Here B a . a (Rj, r„) are the parameters of anisotropic exchange, aj = Xj,yj,Zj are the local 
principal axes for the j-th impurity "spin" and a = x,y,z are the principal axes of magnetic 
matrix. s(Rj) = {s Xj (Rj), s Vj (Rj), s Zj (Rj)} is the effective impurity spin located at Rj and 
S(Rj + r„) = {S x {Hj + r„), Sy(Hj + r„), S z (Rj + rv)} is the spin of the host matrix located 
in the vicinity of j-th impurity. 

We can express impurity spin components in terms of creation and annihilation operators: 

s Zj (Rj) = c+cj - 1/2, s+(R i ) = c+, s_(Rj) = c j: 

s±(Rj) = s Xj (Rj) ± is yj (Rj) (A2) 

with the following commutation relations: 

(cjf = (cjf = °. c J c t + c l c i = 
[ C j, c k ] = 0, [cj, 4} = (1 - 2c^Cj)6 jk . 

In the case of only coherent motion of the host spins S(Rj + r„) = (VoM s //ry)m, 
where Vo is the volume of one host spin. Substituting (|A2|) into ( |Al| ) and writing 
W ex = (H ex )/V, we obtain Eq.flD with w j;a = (V /ti 2r yM s ) J2„ B z a (Rj, r v ), e j>a = 



(V /2h 2 jM s ) B x . >a (Rj, t v ) - iB y JRj. r p 



, a = x,y,z. 
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Captions: 
Fig.l 

Illustration of single domain grain and one impurity ion. 
Fig.2 

Field component of the magnetization versus scaled time. 



9 



impurity 

level 

splitting 




1.0-i 




time, yH K t 



